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Summary 

Let S be the multiplicative semigroup of q x q matrices with positive entries such that every 
row and every column contains a strictly positive element. Denote by (X n ) n >i a sequence of 
independent identically distributed random variables in S and by X^ n ' = X n ■ ■ ■ X\, n > 1, the 
associated left random walk on S. We assume that (X n ) n >i verifies the contraction property 



P([J[X 



(«) 



G S°] ) > 0, 

n>l 

where S° is the subset of all matrices which have strictly positive entries. We state conditions on 
the distribution of the random matrix X\ which ensure that the logarithms of the entries, of the 
norm, and of the spectral radius of the products X^ n \ n > 1, are in the domain of attraction of a 
stable law. 
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I. STATEMENT OF THE RESULT 



Let S be the multiplicative semigroup of q x q matrices with real non negative entries such that 
every row and every column contains a strictly positive element. The subset of S composed of 
matrices with strictly positive entries is a subsemigroup of S denoted by S°. 

Let (ei)i = i t ,_, t q be the canonical basis of the linear space M q . Then a q x q matrix is identified with 
an endomorphism of M q . We denote by (-,-) the canonical scalar product on M q , and we define 
the cones C and C by 

C = {x : x £ lR q ,Vi = 1,... ,q, (x,^) > 0}, C = {x : x G M q , Vi = 1, . . . , q, (x,e;) > 0}. 

If g is a q x q matrix, g* will stand for its adjoint. We have g £ S [resp. g £ S°] if and only if 
g(C) C C and g*{C) C C [resp. either g(C\{0}) C C or g*(C\{0}) C C). 

The product of g and g' in S is denoted by gg', and for x £ C, gx is the image of x under g. Finally 
M q is endowed with the norm || • || defined by 

q 

x £ M q , \\x\\ = ^|(x,ej)|. 

i=i 

Let (X n ) n >i be a sequence of independent identically distributed (i.i.d) random variables (r.v) in 
S defined on a probability space (£1,J-,1P). We consider the left random walk 

X (n \ n > 1, = X( n+1 )(cj) = X n+ i(a;)X( n )(a;). 

Our basic assumption is that (X n ) n >i verifies the contraction property 

(C) iP( (J [I (n) G 5°]) > 0. 

n>l 

The subsemigroup S° is in fact an ideal of S, that is : if g £ S° and </ G S, then g'^ and 55' G S°. 
Consequently S° is stochastically closed for the random walk (X^) n >i. We set 

T{uj) = inf{n : n > 1, G 5°}. 

It is easily shown, Lemma II. 1, that : (C) JP[T < +00] = 1 iP(U n >i[X( n ) G S 10 ]) = 1. 

Our aim is to present conditions on X\ ensuring the distributional convergence to a stable law for 
the sequences of real random variables 

K l [T<n] ln(y,X^x)) n>1 , x,y£C\{0}. 

Denoting by 1 the vector in M q whose all entries equal 1, we point out that the scalar products 
(y,X^x), x,y £ C\{0}, include : 

- the matrix entries : (ej, X^ej), i,j = l,...,q, 

- the norm of the image under X^ of any x £ C\{0} : \\X^x\\ = (l,X^x), 

- the norm |||X( n )||| = (1,X^1} oiX^. 

Closely related to these quantities is the spectral radius A n of the matrix X^ n \ Actually the 
Perron- Frobenius Theorem yields A n > 0, and we shall see that the above mentioned distributional 
convergences also concern the sequence (lnA n ) n>1 . 
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To state our result, one needs the two following real r.v : 

<? <? 
Ni = |||Xi||| = V (ei,Xiej), and Vi = min V(ej,Xiej). 

iVi takes in account the size of the matrix Xi while Vi measures the smallness of its lines. 

Theorem I. Assume that (C) holds and that there exist a real number a, < a < 2, a slowly 
varying function L which is unbounded in case a = 2, and finally some positive constants c + and 
c_ with c + + c- > such that 

(i) lim — — iP[Ai > e u ] = c + , lim — — iP[Ai < e~ u ] = c_, 

u-++oo L[u) u^+oo L{U) 

(ii) lim sup— — JP[Vi < e~ u ] < +oo. 

Then there exist a sequence (a n ) n >i in IR* + with lim n a n = +oo and a sequence (6 n )n>i i n JR such 
that, for any sequences (x n ) n >i and (y n ) n >i of unit vectors of C , the random sequences 

(— (l[T<n]Hyn,X^ n) x n ) -b n ) ) and (— (hxAn-bS) 

converge in distribution to a stable law of index a. 

Observe that Hypothesis (i) means that the real r.v. In N\ belongs to the domain of attraction of a 
stable law of index a, < a < 2, the standard Gaussian case being here excluded since L is assumed 
to be unbounded when a = 2. As it will be seen later on, the hypotheses of Theorem I imply that 
the above considered sequences of random variables have the same distributional behaviour that a 
sum of i.i.d random variables (See § IV). However it is worth noticing that this is not true when 
a = 2 and L is bounded. In fact, to complete the Gaussian case, recall it is proved in [10] that, if 
Conditions (i) and (ii) are replaced by the moment condition 1E[\ In Ai| 2 ] In Vi| 2 ] < +oo, then 
the random sequences of Theorem I converge to a normal law. The method used in [10] is based 
on martingale techniques, and the central limit theorem proved there is also valid when (X n ) n is 
supposed stationary and satisfies suitable mixing conditions. By the way, recall that, in some cases, 
the unnormalized random products (X^) n >i converge in distribution, see [10] [15] [16]. 

Consider the case q = 1. Then Theorem I corresponds to the well-known statement of convergence 
to stable laws for i.i.d random variables : we have S = M*^_, Condition (C) holds, and Condition 
(i) states that InXi is in the above described domain of attraction, (ii) is a consequence of (i). So 
Theorem I gives the expected conclusion. 

The proof of Theorem I is based on the spectral method that was introduced by Nagaev [17], [18] 
and later developped by several authors, see [11]. Although this method has been essentially used 
to prove Central Limit Theorems and their refinements, we mention that Nagaev himself [17] has 
considered the convergence to stable laws, and that his method has been extended to the context 
of dynamical systems, See e.g [8] [2] [3] [4] [9]. 

Section II summarizes some statements of [10], based on Condition (C) and related to the pro- 
jective action defined by g ■ x = for g £ S and unit vector x in C. In Section III, denoting 
by Yfc the adjoint matrix of X^, and setting £(g,x) = ln||^a;||, we show that the distributional 
convergences of Theorem I are valid if, for any unit vector y in C, the same holds for the random 
variables (Yfc_i . . . ii) • y). Since these r.v may be seen as a functional of the Markov chain 
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(Y/i, . . . Y\) ■ y)^, Nagaev's method applies. In fact, it shall be applied to the transition prob- 

ability P of the simpler Markov chain (Y& . . . Y± ■ y)k, and we shall prove, by using the contractivity 
properties stated in Section II, that P satisfies a strong ergodicity condition on a certain Lipschitz 
function space, and finally, by applying the perturbation theory, that the Fourier kernels Pt asso- 
ciated to P and £ inherit near t = the spectral properties of P. 

As usual in Nagaev's method, the previous preparation will show that the desired distributional 
convergence is based on the behaviour of the power of the dominating eigenvalue \{t) of Pt- Actu- 
ally, one of the main arguments is Proposition III. 1 which links X(t) with the characteristic function 
of the r.v £ under the stationary distribution of (Yj,, (l/c_i . . . Y\) ■ y)^. So everything shall turn out 
as in the i.i.d case, provided that £ belongs to the already mentioned domain of attraction. We 
shall see in Section IV that this requirement holds under Conditions (i) and (ii). 

The above relation between the dominating perturbed eigenvalue of the Fourier kernels and the 
characteristic function of the functional under invariant distribution has been already exploited in 
[9] [4], and mentioned under a different form in [11] (Lem. IV. 4'). It is worth noticing that such 
a relation holds whenever the spectral method applies, and that it greatly makes easier the use of 
Nagaev's method when dealing with stable laws excluding the standard Gaussian case ; for instance 
it yields a significant simplification of some proofs in [2] [3]. 

II. CONTRACTIVITY 

II. 1. Preliminaries. We set 

B = C n {x : x G M q , \\x\\ = 1}, and B = C n {x : x G M q , \\x\\ = 1}, 

and we define the adjoint random walk (Y^) n >i of {X^) n >\ by 

Y n = X*, Y W = = Yi ■ ■ ■ Y n , n > 1. 

Lemma II. 1. (C) is equivalent to 1P[T < +oo] = 1. Let uj be such that T(uj) < +oo. Then 
(i) forn> T{uj), X^ 1 \uj) G S° , 

(ii) setting, forn> 1, D n (u) = sup{ 1 [T <„]M ln{y,X {n \uj)x) - ln\\Y^ (u)y\\ : x,y G s}, 

have sup D n (io) < +oo, 
n>l 



we 



(Hi) setting x = 1^ and, for n > 1, D n (uj) = In A n (a;) — In ||Y"^(a;)x|| 



, we have sup D n (u) < +oo. 

n>l 



Proof. Suppose (C) holds. Then there exists k G M* such that p = JP[X^ G S°] > 0. The r.v. 

T' = inf{n : n > 1, X nk . . . X (n _ 1)fe+1 G S°} 

has a geometric distribution with parameter p. Since S° is an ideal, we have T < kT', hence 
IP[T < +oo] = 1. The converse implication is obvious. 

Now let any fixed uj G £1 be such that T{oj) < +oo. Assertion (i) follows from the fact S° is an 
ideal. To prove (ii), it suffices to establish that sup D n (u) < +oo. In the following inequalities, 

n>T(u>) 

one considers any fixed integer n such that n > T(oS). For convenience, oj will be omitted in most of 
the next computations. Let a, b be two strictly positive real numbers such that, for i,j = l,...,q, 
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we have a < (ei, ej) < b, and let x and y be any elements of B. Using ||x|| = ||y|| = 1, we 
obtain for i = 1, . . . , q 

a < (ei,X {T) x} < b. 

(That is, al < X^x < bl for the coordinatewise order relation on M q .) Moreover, using the 
formula (y, X^x) = (lr+i ■ ■ ■ Y n y, X^x), one gets successively 

a (Y T+1 . . . Y n y, 1) < (y, X^x) < b (Y T+1 ■ ■ ■ Y n y, 1), 

| ln{y,X {n) x) -ln\\Y T+1 ---Y n y\\\ < max{| lno|, | ln6|}. 
In particular, with x = 1, this gives 

| ln||y( n )y|| -ln||y T+ i---y n y|| | < max{| lna|, | lnfe|}. These 
two inequalities imply sup D n (oj) < +oo. 

n>T(u>) 

To prove (iii), again consider lo £ Q such that T(cj) < +oo, and recall that, from the Perron- 
Frobenius Theorem, there exists R n (u) G B such that X^ a \u))R n (uj) = K n {uj)R n {uj). With x = R n 
and y = x = hi, Assertion (ii) yields 



l [T < n] \n-{l,X^R n ) - ln\\Y^ X \ 



<D n . 



q 

From (l,X^R n ) = A n \\R n \\ = A n , it follows that 

D n (u) < D n (u) + l [T>n] (uj) | lnA n (w)| +l [T < n] (u)lnq. 
This proves assertion (iii). □ 

We deduce from the above lemma that, for any sequence (a n ) n >i in ]R* + such that lim n a n = +oo, 

we have lim — D n = and lim — D n = a.s. 

n a n n a n 

Consequently, the conclusion in Theorem I for (lpxn] ^ n (lJn, X^x n )) n and (lnA n ) n will hold if the 
same is valid for (In ||y( n )y n ||) n for any sequence (y n )n>i of vectors of B. 

II. 2. Projective action of positive matrices. It is well known that the projective action of 
matrices plays a key part in the study of the asymptotic behaviour of random invertible matrix 
products, cf. [6] for example. As shown in [10], this is also true in the case of positive matrices. 
About this action, we now recall the facts that we shall use throughout, refering to [10] for more 
details and for the proofs. 

Consider the subset C of the ^-dimensional projective space associated with the cone C. In other 
words, C is the set of lines through and some point in C\{0}. These may be represented by 
points of the closed polygon B. An element g G S maps a line in C onto a line in C, and this 
defines its projective action on C. As C is represented by B, the projective action of g moves to 
the action on B defined by 

gx 

g ■ x = j — 77. 

\\9x\\ 

(Recall that gx is the image of x under the linear action of g.) The projective action has the 
following basic properties : if e stands for the identity matrix and g,g' £ S, x £ B, we have 
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e-x = x, (gg') ■ x = g ■ (g' ■ x). 



It is well known [5] that, when B is equipped with the Hilbert distance da, the elements of S have 
a contractive action on B, and that this contractive action is strict for elements of S°. However, 
because the Hilbert distance is unbounded and only defined on B, it is more convenient for our 
purposes to use a bounded distance d on B which have similar properties. This distance, already 
used in [10], is defined as follows. For x = (x\, . . . , x q ) and y = (yi, . . . , y q ) in B, we write 

m(x, y) = sup{A : A G JR + , Mi = 1, . . . , q, Xyi < Xi} = min{y~ 1 Xj : i = 1, ... ,q, yi > 0}. 

1 — s — 

Besides let if be the one-to-one function on [0, 1] defined by ip(s) = . Then, if x,y G B, one 

1 ~\~ s 

has Yh=i x i = J2t=i Vi = 1) thus < m(x, y) < 1, so one may define 

d(x,y) = ip(m(x,y)m(y,x)). 

Proposition II. 1. (cf. [10], § 10) The map d defines a distance on B having the following properties 
(i) sup{d(x, y) : x, y G B} = 1 (ii) if x,y G B, \\x - y\\ < 2d(x, y) 

(iii) the topology of (B,d) is the topology induced on B by the standart topology of M q . 
Moreover, for g G S, there exists c(g) such that 

(iv) ifx,y G B, d(g-x,g-y) < c(g)d(x,y) < c(g) (v) c(g) < I, and c(g) < 1 if and only if g G S° 
(vi) ifg> G S, c(gg') < c{g)c{g') (vii) c(g*) = c{g). 

For any x G B\B and any y G B, we have m(x, y) = 0, so that d(x, y) = 1. Thus 

b\b = u xe B\ B {y :yeB, d(x, y) < 1/2}. 

is an open subset of (B,d). It follows that the topology of (B,d) and the topology induced by 
M q on B do not coincide ; from (ii) the former is finer than the latter. In the sequel, unless 
otherwise stated, when we appeal to topological properties of B and B, we shall assume that these 
sets are endowed with the topologies induced by M q ; the distance d will be only used to express 
contractivity. 

II. 3. Stochastic contractivity. Denote by /i the probability distribution of Y\ = and by /i*™) 
the distribution of = Y± . . . Y n , n > 1. For n > 1, we set 

c(^) = sup{ f s ^^'^ ^{9) :y,y'eB,y^ y'}. 

Since c(-) < 1, we have c(// n )) < 1. Furthermore, the sequence (c(fi^)) n>1 is clearly submulti- 
plicative, so we can define 

k(h) = limc(M (n) )" = inf c{^)^ . 

n n >l 

Using Assertion (v) in Proposition II. 1, it is easily shown that (C) is equivalent to < 1. 

Theorem II. 1. Under Condition (C), there exists a r.v. Z\ taking values in B such that (Y^ n ^ •x) J1 
converges a.s to Z\, the convergence being uniform for x G B. The probability distribution v of 
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Z\ verifies v(B) = 1. It is the unique [i-invariant probability distribution on B, i.e. the unique 
probability distribution on B such that, for any bounded continuous function f on B, we have 

j_(j G f{9-x)dii{gj)dv{x) = f_f{x)du{x). 

Proof. Using the contractivity properties of c(-), we see that the sequence of positive r.v. (c(Y^ n •*))„; 
decreases and hence converges almost surely. Under (C), there exists an integer b G IN* such that 
E[c(Y^)} < 1. The independence then yields limsup fc E[c(Y {kb ))\ < lim k (lE[c(Y^}) k = 0. It 
follows from these two facts that lim n c(y( n )) = a.s. Notice that, by means of the subadditive 

ergodic theorem, we can get, more precisely, lim(c(y^ n ^)) ~= k a.s. 

Set fti = {uj : \\m n c{Y^ n \uj)) = 0}. Let uj G For n > T{uj), the polygons K n (u) = Y^(uj)-(B) 
form a decreasing sequence of compact subsets of B, so that K(uj) = n n >±K n (uj) / 0. Moreover, 
for the distance d, the diameter A (a;) of K(ui) is equal to 0. Indeed we have for n > T(oj), 

A(w) < A n (w) = sup{d(Y^ (uj) • X, Y^ n \uj) -y):x,y£B}< c(Y^(lu)). 

Define Z x (uj) by setting K(to) = {Zi(u)}. Then Z x G K n {uo) implies d(Y { - n \u) ■ x,Zi(lo)) < 
c(Y( n \uj)), and (ii) in Proposition II. 1 yields the desired convergence. 

Now denote by v the law of Z\. Set Z2 = lim n (Y2 • • • Y n ) ■ x a.s. Clearly Z2 has the distribution 
v, and we have Y± ■ Z2 = Z\ a.s. This gives the ^-invariance of v. Let 1/ be any /x-invariant 
distribution on B. Then, for any continuous bounded function / on B and n > 1, we have 

J_IE[f(Y^ ■ x)]du'(x) = J_f(x)du'(x). Thus TE\f(Z x )\ = f(x)di/ (x). Hence v = v' . □ 
III. FOURIER KERNELS 

III.l. Definition and link with our distributional problem. Recall that our aim is to 
study the distributional behaviour of the sequences (In ||Yi • • • V^y||) n>1 , y G B (cf. the end of § 
II.l). However, since ((Y n . . . Y\ y)) n> i is a Markov chain and Y n . . . Y\y has the same distribution 
as Y\ ■ ■ ■ Y n y, it is more convenient to consider (In \\Y n . . . liy||) n>1 . So we introduce the new left 
random walk on S, 

Y^=Y n ...Y u n>l, y(°)=e. 

For yo £ B, consider the sequence of r.v. in B defined by (y( n ) • yo) n >o- ^ ^ s easn y checked that it 
is a Markov chain on B starting at yo and associated with the transition probability P defined by 

Pf(x)= I f(g-x)d^g), 
Js 

where x G B and / is a bounded measurable function on B. Theorem II.l shows that v is the 
unique P-invariant distribution. Finally, for g G S and x G B, define 

£(g,x) = ln||sfx||. 

The function £ is connected with the projective action of S on B by the additive cocycle property 

t(gg', x) = £{g, g' ■ x) + x) (g, g' G S, x G B). 
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This property shows that, for any y £ B and n > 1, we have 

n 

(*) in||Y-w y || = = Y.^y {k ~ l) -y)- 

k=l 

With the function £ and the transition probability P, we associate the Fourier kernels Pt, t € M, 



x£B, P t f(x)= { e it ^f(g-x)d f i(g), 
Js 



is 

with / as above. The Markov property implies that for n > 1, y £ B and t £ M (see e.g [11]) 

(**) lE[e mn ^ (n) ^} = P t n l(y), where 1 = %. 

This basic relation shows that limit theorems for the sequence (^(Y^ n \ x)) n>l may be deduced from 
the asymptotic behaviour of the iterates of the operators Pt acting on a suitable Banach space. 
This is the main idea of the spectral method. In Sections III. 2-4 below, we shall prove that P 
satisfies a strong ergodicity property on the usual space of Lipschitz functions on B, and we shall 
apply the standard operator perturbation theorem to the Fourier kernels. 

III. 2. A strong ergodicity property for P. We denote by L the space of all complex-valued 
functions / on B such that 



m 



»(/) = S M J {x x >) ■■x>x'£B,x^ x'} < +oo. 
Since the distance d is bounded, the elements of C are bounded, so we can equip C with the norm 
/£A ll/L = ll/ll« + m(/), with ||/||„ = sup{|/(aO|:a;€B}. 

Then (£, \\ ■ \\ c ) is a Banach space. Notice that the functions in C may be discontinuous on B w.r.t 
the induced topology of M q , see the remark following Proposition ILL We still denote by || • || £ 
the operator norm on C, and II stands for the rank one projection on C defined by : 11/ = v(f)l- 

Theorem III.l. Under Condition (C), for any kq 1[ ; there exists C > such that, for all 

n > I, we have \\P n - U\\ c < C/eg. 

Proof. We follow [11]. For x, x' G M, x / x', we have 



d(x,x') J d(g-x,g-x') d(x,x') 

so P n f G C and m(P n f) < m(f) c{^). Since ||P/||„ < ||/||„, P acts continuously on C. 
Now set H = Ker(z^) n C. Since v is P-invariant and defines a continuous linear functional on 
C, H is a closed P-invariant subspace in £. Moreover, when restricted to H, the semi-norm m is 
equivalent to the norm || • || £ : more precisely, if h £ H, we have 

m{h) < \\h\\ c < (2sup{d(y,y , ) : y , y' £ B} + l)m(h) <3m(h), 
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the second inequality being deduced from the fact that, if v(f) = 0, there exist x\,X2 G B such 
that Re /(xi) = Im f{x 2 ) = 0. Let / G £. Since / - n(/) G F, we have P n (/ - n(/)) G # for all 
n > 1. Hence 

||P"(/ - n(/))|| £ < 3m(r(/ - H(/))) < 3c(^)m(f - = 3 C (^ )m(/ ) < 30(^)11/^. 

Finally, under (C), we have lim n c(/x^) « = k(^) (§ II. 3). This gives the desired statement. □ 

III. 3. The Fourier kernels near 0. To apply the perturbation theory near t = to the Fourier 
kernels Pt, we have to show that Pt is a bounded operator of C, and to study \\Pt — P\\ c when 
t — > 0. For that, we shall need the following notations. For g G S, define 

\\g\\ = sup{||stx|| : x G B}, v(g) = inf{||gar|| : x G B}, and £{g) = | In \\g\\ \ + \ lnv(g)\. 
(Notice that v(g) > 0.) Finally set e(t) = J s mm{\t\£(g), 2}d/j,(g), and observe that lim t ^. £ (*) = 0. 
Theorem III. 2. For t G M, Pt defines a bounded operator of C, and \\P t — P\\ c = 0(e(t) + \t\). 
Proof. Recall Pt is associated to P and £(g,x) = In \\gx\\ (g G S, x G B). 
Lemma III.l. For g G S and z,x,y G B such that d(x,y) < 1, we have 
\t(9,z)\<i(9), \Z(9,x)-t(g,y)\<2hi- ' 



1 - d(x,y) 

Proof. The first inequality is obvious. The second one is Assertion (ii) of Lemma 5.3 in [10], 

for completeness we reproduce the proof here. Let x = (xi, . . . ,x q ), y = (yi, . . . , y q ) G B and 

q q Q Q 

g = [gij]i,j=i,..., q . Then \\gx\\ = J^J^gijXj > m(x,y)^^g ij y j = m(x,y)\\gy\\. As d(x,y) < 1, 

i=l j=l i=l j=l 

the number m(x,y) and m(y,x) are in ]0, 1]. Consequently, the symmetry in x and y yields 

||gx|| 1 

m(x,y) < j — rr < — -, and 

\\gy\\ m(y,x) 

\((9,x) - ^{g,y)\ <max{- In m(y,x),- In m(x,y)} < - \nm(y, x) - \nm(x, y) 

, -\ (AI \ \ , l + d(x,y) 
= -]n<p (d(x,y)) = ln^— [( — y 

1 . 1 + t 

in 

1 - 1 l-t 1-t 2 



For t G [0, 1[, 2 In 7 - In : = In f ^ °> thus follows. □ 



Set A t = P t - P. For / G £ and s G £, we have A t /(x) = ^ (e^ (9 ' x) - l)f(g • x)dfi{g). 
Before we proceed, notice the inequality : Vu, v G -R, |e m — e™| < min{|-u — u|, 2}. Thus 



|A t /(x)| < y mm{\t\£(g), 2} |/( 5 • < <t)\\f\\ u . 

So HAt/IU < e(t) ||/|| tt . Now for x,y G g, write = A(a,y) + with 



e it£(g,x) _ e it£{g,y) 

A{x,y) 



d(x,y) 



f(g-x)d^g) and B(x,y) = j(e^) - 1) - ^ ^ ' v) d^{g) . 
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If d(x,y) > 1/2, we have 

\ e itt( 9 ,x) _ e «to,y)| < mindtll^^) -S(g,y)\,2} (2d(x,y)) < 4min{|t| 1} d(x, y), 
while, for d(x,y) < 1/2, the inequality of Lemma III. 1 gives 

\ e it£{g,x) _ e itC(ff,y)i < 21*1 In ^ < 2C|t|d(x,y), 

l-d(x,y) 

withC = sup{iln T ^ I : < u < 1/2} < +oo. From that, we obtain \A(x,y)\ < (4e(t)+2C\t\) 
Otherwise, since c{g) < 1, 

M _ „ ,&^!| ^» < ,»</)*). 

So m(A(/) < (4e(t) + 2C\t\)\\f\\ u + e(t)m(f), therefore ||A t || £ < 4e(t) + 2C\t\. □ 

III. 4. Spectral properties of P t near t = 0. The following perturbation theorem extends the 
spectral conclusion of Theorem III.l to Pt for t near 0. Let reo be chosen as in Theorem III.l. 

Theorem III. 3. We assume that Condition (C) holds. Let re e]reo,l[. There exists an open 
interval I centered at t = such that, for t £ I, P t admits a dominating eigenvalue X(t) £ <C, with 
a corresponding rank-one eigenprojection H(t), satisfying the following properties : 

limA(t) = l, \\U(t) -U\\ c =0(\\P t - P\\ c ) and sup \\P? - X(t) n U(t)\\ c = 0(re"). 
; tei 

Proof. We only sketch the proof, refering to [7] for the details and using standard notations. It 
follows from Theorem III.l that the spectrum cr(P) of P is contained in {1} U D(0, reo). Since 
t i ^ Pt is continuous (Th. III. 2), there exists to > such that, for |t| < to, we have cr(Pt) C 
D(l, U D(0, re), and a{P t ) n D(l, ^) = {A(t)}, where X(t) is a simple eigenvalue of P t with a 
corresponding rank-one eigenprojectioon LI(t) depending continuously on i. Let T be the oriented 
circle C(0, re). Since (z,i) i— > (z — -Pt)" 1 is continuous on the compact set T x [— to^o]) the formula 
P t n - A(t) n n(t) = 2^ J r z n {z - P t )~ 1 dz leads to the last estimate of Theorem. □ 

The next proposition states a simple expansion for the perturbed eigenvalue \{t). 

Proposition III.l. For t £ I, we have X(t) = fx ® v(e u ^) + 0(\\P t - P\\ 2 C ). 

Proof. Since v defines a continuous linear functional on C and \\Pt — P\\ c — >0 when t— >0, the 
rank-one eigenprojection , defined in Theorem III. 3, is such that ^(n(t)l) — ► ^(111) = 1. So 
one may assume that ^(n(t)l) / for any t £ I, with / possibly reduced. For t £ I, set 
v (t) = (^(n^l))- 1 II(t)l. Then we have X(t)v(t) = P t v(t) and v(v(t)) = 1, therefore 

X(t) = u(P t v(t)) = v{P t l) + v{P t {v{t) - 1)) = M ® i/(e^) + i/((P t - P)(«(t) - 1)), 

the last equality following from v(P(v(t) — 1)) = v{v{t) — 1) = since v is P-invariant. We conclude 
by observing that \\v(t) - 1|| £ = ||v(t) - u(0)|| £ = 0(||II(t) - n|| £ ) = 0(\\P t - P\\ c ). □ 

IV. PROOF OF THEOREM I 
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Let us point out how the results of the previous sections will be used to establish Theorem I. We 
have to study the distributional behaviour of 

n 

\n\\Y^yn\\=J2^ Y ^ {k ~ 1) -yn) 
_ k=l 

for any sequence (y n )n>i of vectors of B. Let y G B, from Theorem II. 1 and the independence of 
Y k and Y^ k ~ l \ the sequence (£(!&, Y^ k ^ ■ y))k converges in distribution to £(Yi, Z2), where Z2 is 
independent of Y\ and has distribution v. From this we may guess that (In || y( n )y n ||) n>0 has the 
same asymptotical behaviour that a sequence of sums of stationary random variables distributed 
as £(li, Z2). This is confirmed by a look at the characteristic functions. In fact, the characteristic 
function of ln||Y( n )y n || is P(t) n l(y n ) whose asymptotic behaviour is, as shown by the spectral 
decomposition of Theorem III. 3, essentially ruled by X(t) n which doesn't depend on y n or on any 
initial distribution. Notice that this can be used as will be done in the sequel to deduce a limit 
theorem from the expansion of \(t) at 0, but also conversely to get an expansion of A(i) at from 
a known limit theorem, see [12] [13]. Now observe that the characteristic function of £(Y\,Z2) is 
H (8) i/(e'*^), which is precisely the first term in the expansion of X(t) in Proposition III.l. So we see 
that if, for a sequence (a n ) n of positive real numbers, we have 

(★★★) \\P_t_ - P\\ 2 = o(-), 

a„ L U 



n>0 



then (p (g) v(ean^)) n is the principal part of the expansion of (A(^)) n , so that In \\Y^y n 

will behave as J2k=i "fc)n>i> wnere (pk)k>i ^ s a sequence of independent random variables 
distributed as £(Y\,Z2). 

Actually we shall show that, under the hypotheses of Theorem I, the law of £(Y±,Z2) is in the 
domain of attraction of a stable law and that (* * *) is verified with the corresponding scaling 
sequence (a n ) n , these two facts lead to the claimed result. Finally observe that Condition (**★) 
is not fulfilled with a n = y/n in the standard gaussian case since it is known that in this case the 
variance of the limit law is not a 2 (£ (Yi, Z2)) [11]. 

Proposition IV. 1. Suppose that Conditions (C) and (i)-(ii) hold. Let F be the distribution 
function of £(Y±, Z2) : 

F(u) = fj,® v{(g,x) : £(g,x) < u} (u £ M). 
Then there exist positive functions p + and p_ defined on JR* + such that, for u > 0, 

. p + (u)L(u) N p-(u)L(u) . , , . , . . 

1 - F(u) = H+y ' y F(-u) = H y ' v ; , with lim p + (u) = c + and lim p(u) = c_. 

Proof. Since £(g,x) = ln||^a;|| and p is the law of Y±, we have F(u) = _P[\\Yix\\ < e u ]du(x). 

J B 

g 

For x G B, we set Nf = \\Yix\\ = < ei,X\ej > X{. As v{B) = 1, one gets for u > 

1-F(u)= f 1P[N? > e u \du{x) and F(-u) = f JP[Nf < e - u ]du{x). 
Jb Jb 

To proceed, we have to compare the tails of Nf, for x £ B, with that of N±. Let uo be such that, 
for u > uo, we have L(u) > 0. For u > uq, and for x 6 B, we set 
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c + {u) = y^lPWi > e u ], c_{u) = y^lPWi < e~ u ], 
L(u) M n J 

„.a a 

L(u) M n J 
Setting m(x) = min Xj, we have the following obvious inequalities m(x)Ni < Nf < Ni, and 

i=l,...,q 

iP[m(x)iVi > e u ] < P[Nf > e u ] < iP[Aq > e u ] 

JP[N! < e~ u ] < P[Nf < e~ u ] < iP[m(x)iVi < e~ u ]. 
Now let s, < e < 1. Suppose that u > is such that e~ £U < m(x), we get 

. .L((l + e)u) i s 
-c + ((l + e)u) „Y ■ /J < c+(x,u) < c+(u), 



L(u) +vv y y ((l + e)u) 
c_H < c-(x,«) < ^-(U " ((1 _ £)uT • 
Since by hypothesis, c + {v) — ► c + and c_(v) — ► c_ when u — > +oo, it follows that 

< liminf c+(x, u) < limsupc + (x,u) < c+, 



(1 + e) a ~ u^+oo 

c_ < liminf c-(x,u) < limsupc_(x, u) < 



«^+oo «^+oo (l-e) a 

Thus lim c + (x,u) = c+ and lim c_(x, u) = c_. 

Lastly JVf < TVi yields iP[iVf > e u ] < iP^ > e M ], while iVf = ||lix|| = £j=i ||Yiej||a;j > Pi gives 
P[Nf < e~ u ] < IP[Vi < e~ u \. Therefore, for any x € B and -u > 0, we have 

u a 

c+(x,u) < c+(u) and c_(x,u) < — -— JP[Vi < e 

L(u) 

and by (i)-(ii), the functions of the variable u in each right term of these inequalities are bounded 
on M + . Now one may conclude. We have 

(1 — F(u)) = / c + {x,u)dv(x) and F(—u) = / c_(x, u)dv(x), 

L(u) Jb L(u) Jb 

and Lebesgue's Theorem implies that these integrals converge to c + and c_ respectively as«^ +oo. □ 

Proposition IV.l means that £(Yi,Z2) belongs to the domain of attraction of a stable law with 
order a, < a < 2, the standard Gaussian case being excluded since, for a = 2, L is supposed to 
be unbounded. 

Let (Hfc)^ be an independent sequence of real r.v distributed as £(Yi, Z2). From Proposition IV.l, 
there exist sequences (a n ) n in ]R* + and (6 n ) n in iR such that { " 1+ '"\^ n ~ bn )n converges in distribu- 
tion to a stable law of order a, see [14]. It is known that lim n a n = +00 and that (a n ) n >i may be 

chosen such that — L(a n ) = 1. 

Proposition IV.2. Suppose that Conditions (C) and (i)-(ii) hold. For any fixed real t, we have 
||P_t_ — P\\ 2 C = o(y t ) when n — > +00. 
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Proof. The notations \\g\\, v(g) and £(g) below have been introduced for Theorem III. 2. We have 
V\ = v(Y\) : indeed, by definition, V\ = minj = i v .. )(? ||Yiej||, so e; £ 5 implies v(Y\) < Vi, and if 
x € B, one has \\Yix\\ = J2j=i \\Yi e j\\ x j > V±, thus v{Yi) > V\. Besides |||d||| = (1)<?1) is a norm 
for q x ^-matrices, while, for g £ S, the quantity \\g\\ corresponds to the matrix norm associated to 
the norm || • || on M q . Since the two previous norms are equivalent, and N\ = |||Xi|||, ||-Xi|| = ||Yi||, 
there exists a constant C > 1 such that C N\ < \\Yi\\ < CN\. 

Now let < (3 < a. Hypotheses (i)-(ii) show that E[(ln + N^] < +oo and JE[(ln~ Vi) 13 ] < +oo. 
From the previous remarks, we deduce that J s ((g) 13 d/j,(g) < +cc (recall fi is the law of Y{). 
Denote by mp the previous integral. If /? < 1, then we have min{|i|£(g), 2} < 2\t\P 1(g) 13 , so 

that e(t) < 2\tfmp. If /3 > 1, then e(t) < 2|t|/ £(g)dfi(g) < 2\t\m p . Thus Theorem III. 2 gives 

||P t - P|| = 0(|i|^) if < a < 1, and ||P t - P|L = O(t) if 1< a < 2. Finally, using —L(a n ) = 1 

and the fact that L is unbounded in the case a = 2, this easily yields the desired statement. □ 

Proof of Theorem I. Let (y n )n>i be any sequences of vectors in B, and let (f> n [resp. ip n ] be 

the characteristic function of ln||y( "^" 11 [resp. of Sl+ -+^"~ 6 " ]. Let 0(t) = /i <g> i/(e ite ) be 
the characteristic function of ^(Yi,^). Let t e M and n G IV* be such that ^- 6 /, and set 
£„(*) = n(^)l(y n ) (cf. Th. III.3). By Theorems III.2-3, one gets \im n £ n (t) = 1. Furthermore we 
have 

n (t) = e - l ^riE[ e i ^ ln|l?( " ) ^ 11 ] = e- U tp^l(y n ) (by (**) in § III.l) 

On 

= e _it ^ A(— ) n l n (t) +0(n n ) (by Th. III.3) 

a n 

= [0(A) + (I))]n £n(t) + ( K n) (p rop . ni.l, IV.2). 

a n n 

Finally, since */> n (i) = e _i *^ <M^) n , one gets 4> n (t) = ip n (t) [1 + o(±)] n l n (t) + 0(n n ), therefore 
lim n 4> n (t) = lim„ ip n (t) . Since ( " 1+ '" + a ^ n ~ bn )n converges in distribution to a stable law of order a, 
the same holds for ( ln||r( "^" ll ~ 6 " ) n . □ 

v tin. ' 
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